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The effect of fractal clusters on magnetic and transport properties of percolative superconductors is 
considered. The superconductor contains percolative superconducting cluster, carrying a transport 
current, and clusters of a normal phase, which act as pinning centers. A prototype of such a 
structure is a high-temperature superconducting (HTS) wire. The superconducting core of first 
generation HTS wires is the conglomeration of superconducting micro-crystallites containing normal- 
phase clusters inside. There are the clusters of columnar defects in superconducting layer, which 
is the current-carrying element of second generation HTS wire. It is found that normal-phase 
clusters can have essential fractal features that affect the vortex dynamics. The fractal dimension 
of the boundary of normal-phase clusters in YBCO films is estimated and the cluster statistics is 
analyzed. Depinning and transport of vortices in fractal superconducting structures are investigated. 
Transition of the superconductor into a resistive state corresponds to the percolation transition from 
a pinned vortex state to a resistive state when the vortices are free to move. It is revealed that a 
mixed state of the vortex glass type is realized in the superconducting system involved. The fractal 
distribution of critical currents is derived and its peculiarities are studied. It is found that there is the 
range of fractal dimension where this distribution has anomalous statistical properties, specifically, 
its dispersion becomes infinite. The current-voltage characteristics of superconductors containing 
fractal clusters are obtained. Dependencies of the free vortex density on the fractal dimension as 
well as the resistance on the transport current are studied. It is found that the fractality of the 
cluster boundary intensifies pinning and thereby raises the critical current. This feature enables the 
current-carrying capability of a superconductor to be enhanced without changing of its chemical 
composition. 

PACS numbers: 74.81.-g; 74.25. Fy; 74.81. Bd 



I. INTRODUCTION 

An essential feature of the clusters of defects in su- 
perconductor lies in their capability to trap a magnetic 
fj UX( i 4 2 4 3 Holding in place the vortices driven by the 
Lorentz force, such clusters can act as effective pinning 
centers.^ 5 This feature can be used in making new com- 
posite superconducting materials of enhanced current- 
carrying capability? '2*2, The cluster structure affects the 
vortex dynamics in superconductors, especially when 
clusters have fractal boundaryAiSiiiii^ii^ In the present 
paper the magnetic and transport properties of com- 
posite superconductors with fractal clusters of a normal 
phase will be considered as well as the phenomena limit- 
ing the current-carrying capability of such superconduc- 
tors will be discussed. 

A prototype of superconductor containing inclusion of 
a normal phase is a superconducting wire. The first gen- 
eration high-temperature superconducting (HTS) wires 
are fabricated following the powder-in-tube technique 
(PIT). The metal tube is being filled with HTS powder, 
then the thermal and deformation treatment is being car- 
ried out. The resulting product is the wire consisting of 
one or more superconducting cores sheathed by a nor- 
mal metal (Fig.^). The sheath endows the wire with the 
necessary mechanical (flexibility, folding strength) and 
electrical (the possibility to release an excessive power 
when the superconductivity will be suddenly lost) prop- 



erties. At present, the best results are obtained for the 
silver-sheathed bismuth-based composites, which are of 
practical interest for energy transport and storage. In 
view of the PIT peculiarity the first generation HTS wire 
has a highly inhomogeneous structured Superconduct- 
ing core represents a dense conglomeration of BSCCO 
micro-crystallites containing normal-phase inclusions in- 
side (see inset in Fig. ^ . These inclusions primarily con- 
sist of a normal metal (silver) as well as the fragments of 
different chemical composition, grain boundaries, micro- 
cracks, and the domains of the reduced superconducting 
order parameter li^iiLii The volume content of a normal 
phase in the core is far below the percolation threshold, 
so there is a percolative superconducting cluster that car- 
ries the transport current. 

The second generation HTS wires (coated conductors) 
have multi-layered film structure consisting of the metal 
substrate (nickel-tungsten alloy), the buffer oxide sub- 
layer, HTS-layer (YBCO), and the protective cladding 
made from a noble metal (silver). Superconducting layer, 
which carries the transport current, has the texture pre- 
set by the structure of an oxide sub-layer. In the super- 
conducting layer there are clusters of columnar defects 
(see Fig. |2J) that can be created during the film growth 
process as well as by the heavy ion bombardmenti^i^ 
Such defects are similar in topology to the vortices, there- 
fore they suppress the flux creep so it is possible to get 
the critical current up to the depairing valueii^ 
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FIG. 1: Schematic representation for the first generation su- 
perconducting wire. The superconducting core is sheathed 
by the normal metal. The transport current creates the mag- 
netic flux concentrated along irregular-shaped rings passing 
through the normal-phase inclusions. The inset shows the 
cross-section of the wire with normal-phase inclusions as well 
as one cluster defined as a set of pinning centers, which are 
united by the common trapped flux. 



II. MAGNETIC FLUX TRAPPING IN 
PERCOLATIVE SUPERCONDUCTORS 

Let us consider a superconductor containing inclusions 
of a normal phase, which are out of contact with one an- 
other. These inclusions may be formed by the fragments 
of different chemical composition, as well as by the do- 
mains of reduced superconducting order parameter. The 
columnar defects of such a kind can readily be created in 
superconducting film at the sites of defects on the bound- 
ary with the substrate (see Fig. |5J ■ A passage of electric 
current through a superconductor is linked with the vor- 
tex dynamics because the vortices are subjected to the 
Lorentz force created by the current. In its turn, the 
motion of the magnetic flux transferred by vortices in- 
duces an electric field that leads to the energy losses. 
In HTS's the vortex motion is of special importance be- 
cause of large thermal fluctuations and small pinning 
energies^, Here we will consider the simplified model of 
one-dimensional line pinning when a vortex filament is 
trapped by the set of pinning centers:- 1 A superconduc- 
tor containing isolated clusters of a normal phase allows 
for effective pinning, because the vortices cannot leave 
them without crossing the surrounding superconducting 
space. The clusters present the sets of normal-phase in- 
clusions, which are united by the common trapped flux 
and surrounded by the superconducting phase. The flux 
can be created both by an external source (e.g., during 
magnetization in field-cooling regime) and by the trans- 
port current (in self-field regime) . The magnetic flux re- 
mains to be trapped in the normal-phase clusters till the 
Lorentz force created by the transport current exceeds 
the pinning force. 

If the relative portion of superconducting phase ex- 
ceeds the percolation threshold, there is a superconduct- 




FIG. 2: Superconducting film with columnar clusters of a 
normal phase. Such a film represents the current-carrying el- 
ement of the second generation superconducting wire (coated 
conductors). 

ing percolation cluster carrying a transport current. In 
such a structure magnetic flux is locked in finite clus- 
ters of a normal phase (see Fig. Efa)). If the transport 
current is passed through the sample, the trapped flux 
remains unchanged as long as the vortices are still held 
in the normal-phase clusters. The larger clusters have to 
allow for a weaker pinning than smaller ones, because the 
larger the cluster size, the more entry points into weak 
links, through which the vortices can pass, are located 
along its boundary. When the current is increased, the 
magnetic flux will break away first from the clusters of 
smaller pinning force, and therefore, of larger size. Mag- 
netic flux trapped into a single cluster is proportional to 
its area A. Therefore, the decrease in the total trapped 
flux A<I> is proportional to the number of the clusters 
of area larger than a preset value A. So that value can 
be expressed with the cumulative probability function 
W (A) = Pr{VA,- < A}, which is equal to the probabil- 
ity to find the cluster of area Aj smaller than an upper 
bound of A: 

-^ = 1-W(A) (1) 

The left hand side of this formula is equal to the rel- 
ative decrease in the total trapped flux caused by the 
transport current of the same amplitude as the depinning 
current of the cluster of area A, and the right hand side 
gives the provability to find the cluster of area greater 
than A in the whole population. 

When the current is increased, the magnetic flux starts 
to break away from the clusters of pinning force weaker 
than the Lorentz force created by the transport cur- 
rent. The vortices, driven by this force, must cross the 
surrounding superconducting space, and they will first 
do that through the weak links, connecting the normal- 
phase clusters (see Fig. |3Ja) where weak links are shown 
by the black curves connecting the white normal-phase 
clusters). In such a system depinning has a percola- 
tive character)22i22i2i because unpinned vortices move 
through the randomly generated channels created by 
weak links. In Fig.|3Jb) the equivalent weak- link- network 
circuit representation of such a system is shown. The 
normal-phase clusters are presented by the cells of the 
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FIG. 3: Equivalent weak- link-network circuit representation 
of superconductor with normal-phase clusters, (a) - the per- 
colation chart (superconducting space is shaded, weak links 
are shown by the curves joining together the normal-phase 
clusters); (b) - the electric circuit (the superconducting loops 
are shown by lines, the crosses on them denote weak links). 



net interlaced by the superconducting loops. Each of the 
loops contains weak links (shown by the crosses), which 
join the adjacent cells and so enable the vortices to pass 
from one cluster to another. 

Weak links form readily in HTS characterized by an 
extremely short coherence Iengthi22i25i2& Various struc- 
tural defects, which would simply cause some additional 
scattering at long coherence length, give rise to the weak 
links in HTS. On a mesoscopic scale twin boundaries 
are mainly responsible for weak link existenceiSLSiSLSi 
Twins form especially readily in YBCO superconductors 
inasmuch as their unit cell is only close to the orthorhom- 
bic one. The twins can be spaced up to several nanome- 
ters apart, so even single crystal may have the fine sub- 
structure caused by twins. Magnetic and transport prop- 
erties of HTS depend strongly on the orientation of twin 
planes with respect to the applied magnetic tieldmi^ 
The flux can easily move along the weak links formed by 
twins. 33,34 "* 5 - 36 At last, on a macroscopic scale there arc 
manifold structural defects which can form weak links: 
that may be grain or crystallite boundaries as well as 
barriers arising from the secondary degrading the non- 
stoichiometric crystal into the domains with a high and 
low content of oxygeniSSiiSLiSSiSLii Moreover, a magnetic 
field further reduces a coherence length, thus resulting in 
more easy weak links formation4i In conventional low- 
temperature superconductors weak links can be formed 
due to the proximity effect at the sites of minimum dis- 
tance between the next normal-phase clusters. 

As soon as the transport current is turned on, this one 
is added to all the persistent currents, maintaining the 
constant trapped magnetic flux. Each of these currents is 



circulating through the superconducting loop around the 
normal-phase cluster wherein the corresponding portion 
of the magnetic flux is locked. The loop contains weak 
links that join the adjacent normal-phase clusters. As 
the transport current is increased, there will come a point 
when the overall current flowing through the weak link 
will exceed the critical value, so this link will switch into a 
resistive state. During this process the space distribution 
of the currents throughout the superconducting cluster is 
changed in such a way that the resistive subcircuit will 
be shunted by the superconducting paths where weak 
links are not damaged yet. Magnetic field created by 
this re-distributed transport current acts via the Lorentz 
force on the current circulating around the normal-phase 
cluster. As a result, the magnetic flux trapped therein 
will be forced out through the resistive weak link, which 
has become permeable to the vortices. 

Thus, whatever the microscopic nature of weak links 
may be, they form the channels for vortex transport. 
It appears that according to their configuration each 
normal-phase cluster has its own value of the critical cur- 
rent, which contributes to the overall statistical distribu- 
tion. By the critical current of the cluster we mean the 
current of depinning, namely, such a current at which 
the magnetic flux ceases to be held inside the cluster of 
a normal phase. When a transport current / is grad- 
ually increased, the vortices will break away first from 
clusters of small pinning force, and therefore, of small 
critical current. Thus the decrease in the trapped mag- 
netic flux is proportional to the number of all the normal- 
phase clusters of critical currents less than a preset value. 
Therefore, the relative change in the trapped flux can 
be expressed with the cumulative probability function 
F (I) = Pr {Vij < /} for the distribution of the critical 
currents of clusters: 



= F(I) 



(2) 



The critical current distribution F = F(I) is related 
to the cluster area distribution W = W(A), because the 
cluster of a larger size has more weak links over its bound- 
ary with the surrounding superconducting space, and, 
consequently, the smaller depinning current. 



III. GEOMETRIC PROBABILITY ANALYSIS 
OF THE DISTRIBUTION OF WEAK LINKS 
OVER THE CLUSTER BOUNDARY 

In order to find out the relationship between the dis- 
tribution of the critical currents of the clusters and the 
distribution of their areas, the distribution of entry points 
into weak links over the boundary of a normal-phase clus- 
ter should be analyzed. We will examine the geomet- 
ric properties of the normal-phase clusters in its cross- 
section by the plane where the transport current flows. 
The vortices are moving transversally to this plane so in 
order to leave the normal-phase cluster they have to get 
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over the barrier at the boundary of that cross-section. 
We will also suppose that the boundaries of the cluster 
cross-section are statistically self-similar. 

The problem of exit of a vortex from a normal cluster 
represents the two-dimensional analogue of a problem of 
a random walk particle reaching a border42ii£ At the 
same time, unlike the classic problem on the distribu- 
tion of the exit points, here the boundary of the area 
is not absorbing all over. In order to leave the clus- 
ter the vortex has to enter into one of the weak links, 
which are randomly arranged along the cluster perime- 
ter. There are only discrete absorption points, which are 
located just at the sites where weak links are going out 
on the cluster boundary. In other words, these points 
are the points of the entry of vortices into weak links, 
or simply, entry points. For simplicity it will have been 
supposed that after the vortex reaches the entry point, 
it passes all the way between two adjacent normal-phase 
clusters without being trapped inside the weak link itself. 
Here the magnetic flux is transferred by Joscphson vor- 
tices. The Josephson penetration depth is large enough 
in the considered materials, so the size of the region, 
where the vortex is localized, much exceeds the charac- 
teristic length of all possible structural defects that can 
occur along the transport channel. Thus the probabil- 
ity that such a vortex, driven by the Lorentz force, will 
be trapped in passing through a weak link is very small. 
This assumption agrees well with the results of research 
on the magnetic flux motion along weak links^S^^ in- 
cluding twins At the same time, it allows us 
to highlight the role played by the cluster boundary in 
the magnetic flux dynamics. 

Let us consider the distribution of entry points over 
perimeter of a normal-phase cluster. Generally, this dis- 
tribution varies from one cluster to another, so that each 
normal-phase cluster has the entry point distribution 
function ip (I) of its own, which belongs to some func- 
tion class n. Here I is the co-ordinate measured along 
the cluster perimeter, so I £ (0, P) . In this context the 
functions of class Q are random elements of the statis- 
tical distribution. The probability distribution of func- 
tions ip (I) over all the clusters can be characterized by 
the functional Pr{ip (/)}, which is equal to the probability 
of finding a given function ip (I). 

In the most general way the geometric probability anal- 
ysis of the entry point into weak link distribution can be 
carried out by means of path integral technique^ The 
probability that there is the function ip (I) in class may 
be expressed by the path integral 

Pi{ip(l)\n}= J Dip(l)Pr{ip(l)} 
(fi) 

Therefore, the most probable function of entry point 
distribution is the mean over all functions of class Q 



* (I) = (I) = J Dip (I) iP (I) Pr {iP (I)} (3) 

(O) 



The path integral Fourier transform on the probability 
functional Pr{ip (I)} represents the characteristic func- 
tional 



H[k(l)\ 



J Vip (I) exp (if dlk (I) ip (/)) Pr {ip (/)} 

I Zty(0 Pr{iP(l)} 
(n) 



(4) 

where k — k (I) are the functions of a reciprocal function 
set, and integration in the kernel exp (i § dlk (I) ip (I)) is 
carried out over the cluster perimeter. 

The characteristic functional is the path integration 
analog for the usual moment-generating function. The 
probability functional Pr {ip (I)} can be written as an in- 
verse path integral Fourier transform on the characteris- 
tic functional 



Pr {ip (I)} = J Dk {I) exp (-i j> 



dlk(l)ip(l) H[k{l)} 



where the path integration is carried out on the reciprocal 
function space. 

In the simplest case, when all the clusters are of 
equal entry point distribution, which coincides with the 
most probable one of Eq. J3} , the probability functional 
Pv{ip(l)} is zero for all ip (I) that differ from 'J (I), 
whereas Prj'I'^)} = 1. Therefore, the characteristic 
functional of Eq. 0} becomes 



H [k (I)} = exp (i j> dlk (I) * (I) 



(5) 



If all entry points had fixed co-ordinates lj instead of 
the random ones, their distribution would be ip (I) = 
[3 X)j=i ^ Q ~ h) -i where N is the number of entry points 
along the cluster perimeter, 5 (I) is Dirac delta function. 
The constant (3 is being chosen to normalize the distri- 
bution function ip (I) to unity, so that j3N = 1. 

Now suppose that all the points of entries into weak 
links are randomly distributed with uniform probabil- 
ity over the cluster perimeter, so the probability to find 
any j'-th point within some interval dlj is proportional to 
its length. In that case the characteristic functional of 
Eq. takes the form 

N ( N \ 

§ J] dlj exp i/3 § dlk (I) 8 (I - l ) 



H[k{l)] = 



i r si 

pN 



N 

§ n di 3 

i=i 

N 



fUdlj exp U85>(Ji) 
^fljdl.e^ 



(0) 



Expanding the function e l/3k ^ in a power series at 
N >> 1, and taking into account the condition (3N = 1, 
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we may write 
1 

P 



exp 



dl k (I) 



that, after substitution into Eq. J§}, gives 

H[k(l)]=exp(i^- j dlk(l)j (7) 

The found characteristic functional of Eq. Q has the 
form of Eq. JSJ for the function of uniform distribution 
of entry points 



*(!) = 



1 

P 



(8) 



This means that all the clusters have the same 
uniform distribution of the entry points of Eq. for 
which the probability of finding a weak link at any point 
of the perimeter is independent of its position. 

Let us suppose that concentration of entry points into 
weak links per unit perimeter length n — N/P is constant 
for all clusters, and all the clusters are statistically self- 
similar. In this case the mean number of entry points N 
along the cluster perimeter is proportional to its length: 



(9) 



N = <bn(l)dl = nP 



Next step will consist in finding the relationship be- 
tween the size of a cluster and its critical current. The 
pinning force corresponds to such a current at which the 
vortices start to break away from the cluster. As the 
transport current is increasing, the Lorentz force, which 
expels the magnetic flux, increases as well. The vor- 
tices start to leave the normal-phase cluster when the 
Lorentz force becomes greater than the pinning force. 
At the same time, growing in current will result in 
re-distribution of the magnetic flux, which will pene- 
trate deeper and deeper into a transition layer on that 
side of the surrounding superconducting space where the 
Lorentz force is directed (see Fig.^b)). In order to leave 
the normal-phase cluster, vortices have to reach the en- 
try points into weak links. The exit of the magnetic flux 
can be considered as the result of random walks of vor- 
tices driven by the Lorentz force, which is pushing them 
into weak links. The similar approach has been success- 
fully applied in analyzing the magnetic flux penetration 
in SQUID arrays44 The mean number of the entry points 
N available on the cluster perimeter gives the probabil- 
ity measure of the number of the random walk outcomes, 
which are favorable for the vortex to go out. In the case 
of the uniform entry point distribution, from Eq. it 
follows that N oc P, so the perimeter length also repre- 
sents the probability measure of the amount of favorable 
outcomes for vortex to leave the cluster. The more entry 
points into weak links are accessible for random walk vor- 
tices, the smaller is the Lorentz force required to push the 



flux out. Hence, we may write the following relation be- 
tween the critical current of the cluster and its geometric 
size: 



T 1 1 

1 cx = oc — 

N P 



(10) 



This expression is true for the simplest case of uniform 
distribution of entry points, which is assumed to be the 
same for all clusters. Such a simplification allows us to 
emphasize that in the case being considered the magnetic 
flux is held in the normal-phase cluster by its boundary. 

Thus, to deal with the distribution function of Eq. QJ, 
the relation between perimeter and area of clusters 
should be studied. It might be natural to suppose that 
the perimeter-area relation obeys the well known geomet- 
ric formula: P cx \f~A. However, it would be a very rough 
approximation, because this relation holds for Euclidean 
geometric objects only. As was first found in Ref.0, the 
boundaries of normal-phase clusters can be fractal, i.e. 
can have non-Euclidean features. The fractal nature of 
such clusters exerts an appreciable effect on the magnetic 
flux dynamics in superconductors^^SiSi 



IV. A BRIEF INTRODUCTION TO FRACTAL 
GEOMETRY 

The notion of a fractal as an object of fractional di- 
mension was first introduced by Mandelbrot in 1967 
and since then it has received a lot of applications in var- 
ious domains of sciences This concept is closely 
connected to ideas of scaling and self-similarity. Self- 
similarity is invariance with respect to scaling; in other 
words, an invariance relative to multiplicative changes 
of scale. Whereas a usual periodicity is invariance with 
respect to additive translations, self-similarity is a peri- 
odicity in a logarithmic scale. 

The simplest examples of self-similar objects are Can- 
tor sets^ Such a set has measure zero and, at the same 
time, has so many elements that it is uncountable. The 
Cantor set may be constructed by such a way. Let us 
draw a line segment, and erase its middle third. Then 
eliminate middle third of each remaining part, and so 
on. The resulting set obtained by endless erasing of the 
middle thirds of remaining intervals is called "Cantor 
dust"i££ The Cantor dust forms a self-similar set because 
it is invariant to scaling by a factor s = 3. But this set is 
not only a self-similar one, but it is a fractal as well. Let 
us see what part of the set falls within one s-th part of 
the original set, where s is a scaling factor, and ask what 
fraction of the set falls into that portion? In other words, 
how many subsets, each is similar to the original set, are 
there if the length is subdivided into s parts? This num- 
ber of such subsets is equal to N = 2, so one-half of 
the original set falls into one third of the initial length. 
The fractal dimension is defined as the logarithm of the 
subset number divided by the logarithm of the scaling 
factor: D = bxN/hxs = In 2/ In 3 = 0.631. This formula 
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represents the relation between the subset number and 
the scaling factor: N = s D . Thus the fractal dimension 
of the Cantor dust is less than unity. This fact reflects 
its dust-like consistency compared to a usual line. An 
Euclidean line is such a set that if we change the length 
scale, we recover the same set of points. Hence the fractal 
dimension of a line coincides with the topological dimen- 
sion, and both of them are equal to unity. Generally, a 
fractal object has a fractional dimension. A fractal set is 
such a set for which the fractal dimension strictly exceeds 
its topological dimension^ 

The Cantor dust gives an example of determinate frac- 
tal that is self-similar by deterministic construction. But 
the same approach can be applied to the stochastic frac- 
tals, which are statistically self-similar^ The fundamen- 
tal feature of any fractals, both determinate and the 
stochastic ones, is that its characteristic measures obey 
a scaling law that includes an exponent called the fractal 
dimension. 

Fractal dimension of statistically self-similar curve can 
be estimated by covering the curve by square grid. The 
size of the grid cell L x L should be small enough for 
measuring all the bend of the curve. If the whole curve 
(the boundary of the normal-phase cluster cross-section) 
goes in square of side L max , then each of N(L) self-similar 
subsets of that curve will go in square of side L — L max /s, 
where s is scaling factor. Let us suppose that the curve is 
self-similar at any scaling factor. The minimum number 
of such squares of side L needed to cover the explored 
curve is equal to 

N(L) = s D = (^pY ocL- D 

where D is the fractal dimension of the curve (so-called 
coastline dimension^ So we have the formula that de- 
termines the fractal dimension 



In (L max /L) 

It is obvious that the magnitude of D is always below the 
topological dimension of a plane, which is equal to two. 

The ratio between the size of measuring cell L and 
the big square side L max sets the unit of measurement of 
length. At £ max = 1 the expression for the fractal dimen- 
sion of Eq. JTTJ takes the form D = —ln(N{L))/ln(L). 

Thus, in order to find the fractal dimension it is neces- 
sary to count the total number N(L) of non-empty cells 
of area L 2 , which cover the curve. To attain the highest 
precision of estimate we have to get wide range of L and 
to average obtained values of D over all the configura- 
tions. 

The length of topologically one-dimensional fractal 
curve of size L max , measured with an accuracy of L, is 
equal to a product of the number of cells, which cover 
that curve, and the linear size of the measuring cell 

P = NL= (h^ D L = |frr (12) 



In other words, if the whole curve goes in the scale of 
length i max , then in the smaller scale of L = L max /s (i. 
e. for the yardstick of a smaller size) the same curve will 
consist of N = s D pieces of length L. As is seen from 
the formula of Eq. 1|12(1 . the length of a fractal curve is 
not well defined, because its value depends on the ac- 
curacy of the measurements, and even diverges as the 
yardstick size L is reduced infinitely. At the same time, 
the area restricted by any closed (including fractal) curve 
is well-defined finite quantity, which is proportional to 
its squared linear size: A cx L^ax- According to ex- 
pression of Eq. (|12[1 . the length of the perimeter of the 
figure, formed by the closed fractal curve, is proportional 
to its linear size raised to a power of fractal dimension: 
P cx imax- Hence it follows that the perimeter of fractal 
figure and enclosed area have to obey a scaling law 

P X ' D cx A 1 ' 2 (13) 

The perimeter-area scaling law of Eq. I|13|) expresses the 
generalized Euclid theorem about measures of similar fig- 
ures, stating that the ratios of corresponding measures 
are equal when reduced to the same dimension^ This 
assertion is valid both for Euclidean figures and for the 
fractal ones. 

The fractal approach has been found to be most useful 
in an investigation of inhomogeneous materials. There 
are many possibilities both the determinate fractals and 
the stochastic ones to be formed in composite supercon- 
ductors. As an example of the first kind the multi-layered 
structures prepared by electron-beam deposition of su- 
perconductor (Nb) and normal metal (Cu) layers with 
fractal stacking sequence on sapphire substrates can be 
mentioned^ In order to obtain stochastic fractal clus- 
ters, it is essential that the process like the diffusion- 
limited aggregation would take place in the course of the 
synthesis of material^ The similar process can be real- 
ized, for instance, when thin films are evaporated. So, in 
Ref. |57| the films of fractal structure has been grown by 
vapor deposition of gold on silicon substrate with silicon 
nitride buffer sub-layer. It is worthy of note that porous, 
random, or highly ramified clusters do not necessarily 
all are fractals. A fractal cluster has such a property 
that its characteristic measures (in what follows - the 
perimeter and the enclosed area) have to obey the cer- 
tain scaling law that includes an exponent named fractal 
dimension . 52 ' 55 

The fractal clusters can be also formed in such highly 
inhomogeneous materials as high-temperature supercon- 
ductors. So, the fractal dissipative regime has been ob- 
served in high-resolution measurements of dynamical re- 
sistance of BSCCO and BPSCCO composites containing 
normal-phase inclusions of Ag^i The fractal properties 
of the normal-phase clusters contained in YBCO films, 
which were prepared by magnetron sputtering on sap- 
phire substrate with a cerium oxide buffer sub-layer, have 
been found in Ref. Percolation clusters are another 
example of fractals in superconductorsi^iS Although, 
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mathematically, the percolation cluster is a fractal at 
the threshold point only, the fractal approach works well 
for any clusters which have a scaling featurep^SiiSS In 
that case the normal-phase clusters may be formed by 
inclusions of different chemical composition, as well as 
the domains of reduced superconducting order parame- 
ter can act as such clusters. The existence of the fractal 
inclusions of this kind can be demonstrated by the frac- 
tal dissipation, which has been observed in non-textured 
polycrystalline YBCO and GdBCO bulk samples»±i The 
fractal structure of clusters near the percolation thresh- 
old in epitaxial YBCO films has been fully considered in 
Ref. |6jJ Of special interest are the works of Refs. IToIImI 
where the fractal penetration of magnetic flux in thin 
HTS films has been investigated by magneto-optical tech- 
nique. Epitaxial T1BCO films were grown by magnetron 
sputtering on SrTiOa substrates, and YBCO films were 
prepared by pulsed laser deposition on NdGa03 sub- 
strates. The cluster structure of such films is clearly visi- 
ble on the Atomic Force Microscopy picture published in 
Ref. IToL whereas the magnetic flux, penetrating into the 
sample from the outside, has well-definite fractal front. 

V. FRACTAL GEOMETRY OF 
NORMAL-PHASE CLUSTERS 

The scaling law of Eq. ((13(1 furnishes a clue to find 
relation between the critical current of the cluster and 
its geometric size. Using the formula of Eq. pOjl. which 
links the cluster size and its current of depinning, can 
we get the following expression: / = uA~ D / 2 , where a 
is the form factor, and D is the fractal dimension of the 
cluster boundary. In the general way the cluster area dis- 
tribution can be described by gamma distribution, which 
is appropriate for the most part of practically realizable 
cases 42t&M&t&&i This distribution is characterized by the 
following cumulative probability function: 

W(A) = (T(g + l))- 1 1 (g + l,^J (14) 

where T (y) is Euler gamma function, 7 (u) is the incom- 
plete gamma function, Aq and g are the parameters of 
gamma distribution that control the mean area of the 
cluster A = (g + 1)Aq and its variance a\ — (g + 1) Aq. 

In accordance with starting formulas of Eq. JQ) and 
Eq. J2J, gamma distribution of cluster areas of Eq. 114(1 
gives rise to the critical current distribution of the form: 

F(i) = (r(g + l))- 1 r(g + l,G l - 2 ? D ) (15) 

where G = (0 9 / (9 g+1 - (D/2) exp (0) T (g + 1, 9))) 2/D , 
= g + 1 + D/2, r (y, z) is the complementary incom- 
plete gamma function, i = I / I c is the dimensionless elec- 
tric current, I c — a (AqG) D ' 2 is the critical current of 
the transition into a resistive state. The found cumula- 
tive probability function of Eq. ((15(1 allows to derive the 



probability density f(i) = dF/di for the critical current 
distribution: 

m = rW-X^M exp (-Gi-*l°) (16) 

This distribution allows us to fully describe the effect 
of the transport current on the trapped magnetic flux 
taking into account the fractal properties of the normal- 
phase clusters. Let us note that the probability density 
for the critical current distribution of Eq. 1 (16(1 is equal 
to zero at i = 0, which implies the absence of any contri- 
bution from negative and zero currents. This will allow 
us to avoid any artificial assumption about the existence 
of a vortex liquid, having finite resistance in the absence 
of transport currents due to the presence of free vortices: 
r (i — > 0) 7^ 0. Such an assumption is made, for example, 
in the case of normal distribution of critical currents^ 
In order to clear up how the developed approach can 
be used in practice, the geometric probability analysis of 
electron photomicrographs of superconducting films was 
carried out. For this purpose electron photomicrographs 
of YBCO film prepared by magnetron sputtering on sap- 
phire substrate with a cerium oxide buffer sub-layer have 
been scanned. The normal phase has occupied 20% of 
the total surface only, so the transport current can flow 
through the sufficiently dense percolation superconduct- 
ing cluster. The perimeters and areas of clusters have 
been measured by covering their digitized pictures with 
a square grid of spacing 60x60nm 2 . The results of the 
statistical treatment of these data are presented in Fig. 0] 
The primary sampling has contained 528 normal-phase 
clusters located on the scanned region of a total area of 
200 /1m 2 . The distribution of the cluster areas is fitted 
well to exponential cumulative probability function, as is 
shown on the histogram in the lower inset of Fig.^J The 
number of clusters that fall within the assigned rank is 
plotted on the ordinate of this graph; the rank number is 
plotted on the abscissa. A high skewness (1.765) as well 
as the statistically insignificant (5%) difference between 
the sample mean area of the cluster and the standard 
deviation also attests that there an exponential distribu- 
tion of the cluster areas holds true. This distribution is 
a special case of gamma distribution for which g = 0, so 
the cumulative probability function of Eq. 1(14(1 can be 
simplified to the following form: 

W(A)=l-exp(-^j (17) 

The exponential distribution has only one characteristic 
parameter - the mean cluster area A. 

The obtained data allow us to find the perimeter-area 
relation for the normal-phase clusters. All the points in 
Fig. 0] fall on a straight line in double logarithmic scale 
with correlation coefficient of 0.929. Accordingly to the 
scaling relation of Eq. 1(13(1 . the slope of the regression 
line gives the estimate of fractal dimension of the cluster 
perimeter, which is equal to D — 1.44 ± 0.02. The graph 
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FIG. 4: Perimeter-area relationship for the normal-phase 
clusters with fractal boundary. The chart is built for the 
primary sampling containing 528 normal-phase clusters. The 
insets represent histograms of sampling for corresponding dis- 
tributions. 
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FIG. 5: Crossover between the fractal and Euclidean regimes 
on dependence of the magnetic flux trapped in fractal clusters 
of a normal phase on a transport current. The solid line 
shows the decrease in trapped flux for the fractal clusters of 
boundary dimension D = 1.44; the dotted line corresponds to 
the case of Euclidean clusters (D — 1); step line is the sample 
empirical function of critical current distribution. 



in Fig. 01 shows that the scaling law for perimeter and 
area, which is inherent to fractals, is valid in the range 
of almost three orders of magnitude in cluster area. The 
perimeter-area scaling behavior means that there is no 
characteristic length scale in the range of two orders of 
linear size of the normal-phase cluster. Whatever the 
shape and size of the clusters may be, all the points fall 
closely on the same straight line in logarithmic scale; so 
that there are no apparent kinks or bends on the graph. 
This point that the found value of the fractal dimension 
differs appreciably from unity engages a great attention. 
What this means is that the fractal properties of the clus- 
ter boundary are of prime importance here. 

The geometric probability properties of the normal- 
phase clusters are responsible for main features of the 
critical current statistical distribution. Now, knowing the 
fractal dimension of the cluster boundaries, the change 
in the trapped magnetic flux caused by the transport 
current can be found with aid of Eq. (J2J. The expo- 
nential distribution of the cluster areas of Eq. lfT7|l gives 
rise to the exponential-hyperbolic distribution of critical 
currents 

F(zHexp(-(^^ + V 2 /^ (18) 

which follows from Eq. Ijl5|l at g = 0. The effect of a 
transport current on the trapped magnetic flux is illus- 
trated in Fig.[S]for the case of Euclidean clusters (dotted 
curve), and for the clusters of found fractal dimension 
D = 1.44 (solid curve). 

In order to get the relationship between the dynamics 
of the trapped magnetic flux and the geometric properties 



of the superconducting structure the empirical function 
of the distribution of the critical currents F* — F*(i) 
has been found. First, the empirical distribution func- 
tion W* — W* (A) for the sampling of the areas of the 
normal-phase clusters has been obtained. The value of 
W* (A) was calculated for each order statistic as the rela- 
tive number of clusters of area smaller than a given value 
A. Next, the empirical distribution of the critical cur- 
rents was computed for the same order statistics using 
the formulas: 



F* = 1— W* 
i = {(2 + D) /2) (2+D)/2 (A/A) 




This function, shown in Fig. Elby the step line, gives a 
statistical estimate of the cumulative probability function 
of Eq. i|18[l . As is seen from this figure, both distribu- 
tions coincide well in the range of currents i < 6. Start- 
ing with the value of the current i = 6 the crossover from 
the fractal regime to the Euclidean one is observed. This 
transition into the Euclidean regime is over at large trans- 
port currents, when the magnetic flux changes mainly for 
the breaking of the vortices away from the small clusters 
(as the smaller clusters have the larger pinning force). 
The observed crossover has its origin in the finite reso- 
lution capability of measuring the cluster sizes. When 
estimating the fractal dimension, we have to take into 
account that the resolution of the measurement of any 
geometric sizes is finite. The peculiarity of the topolog- 
ically one-dimensional fractal curve is that its measured 
length depends on the measurement accuracy^ In our 
case such a fractal curve is represented by the boundary 
of the normal-phase cluster. That is why just the statis- 
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tical distribution of the cluster areas, rather than their 
perimeters, is fundamental for finding the critical current 
distribution. The topological dimension of perimeter is 
equal to unity and does not coincide with its fractal di- 
mension, which strictly exceeds the unity. Therefore the 
perimeter length of a fractal cluster is not well defined, 
because its value depends on the yardstick size. On the 
other hand, the topological dimension of the cluster area 
is the same as the fractal one (both are equal to two). 
Thus, the area restricted by the fractal curve is a well- 
defined quantity. 

Taking into account the finite resolution effect, the 
perimeter-area relationship of Eq. l|13fl can be re-written 
in the following form 

P(S)<x8 l - D {A(S)) D/2 (19) 

where S is the yardstick size used to measure this length, 
(1 — D) is the Hausdorff codimension for the Euclidean 
one-dimensional space. This relation holds true when 
the yardstick length is small enough to measure accu- 
rately the boundary of all smallest clusters in sampling. 
When the resolution is deficient, the Euclidean part of 
the perimeter length will dominate the fractal one, so 
there is no way to find the fractal dimension using the 
scaling relation of Eq. I|19|) . It means that if the length 
of a fractal curve was measured too roughly with the 
very large yardstick, its fractal properties could not be 
detected, and therefore such a geometrical object would 
be manifested itself as Euclidean one. It is just the reso- 
lution deficiency occurs at the crossover point in Fig. [S] 
Starting with the cluster area less than some value (which 
corresponds to the currents of i > 6) it is impossible to 
measure all "skerries" and "fjords" on the cluster coast- 
lines, whereas all the clusters of area less than the size of 
the measuring cell (3600 nm 2 that relates to the currents 
of i > 23), exhibit themselves as objects of Euclidean 
boundaries. This resolution deficiency can be also ob- 
served in Fig. 0] some points at its lower left corner are 
arranged discretely with the spacing equal to the limit 
of resolution (60 nm) , because some marks for smallest 
clusters coincide for the finite resolution of the picture 
digitization procedure. 

The fractal dimension was found above by means of 
regression analysis of the whole primary sampling, where 
the very small clusters of sizes lying at the breaking point 
of the resolution limit were also included. In order to 
evaluate how the finite resolution affects the accuracy of 
the estimation of the fractal dimension, all the points, 
for which the resolution deficiency was observed, were 
eliminated from the primary sampling. So the truncated 
sampling has been formed in such a way that only 380 
clusters, for which the resolution deficiency is not ap- 
peared yet, have been selected from the primary sam- 
pling. The regression analysis of the truncated sampling 
is presented in Fig. H3 The least squares treatment of 
perimeter-area data gives the adjusted magnitude of the 
fractal dimension: D = 1.47 ± 0.03. The found value 
virtually does not differ from the previous one within the 
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FIG. 6: Perimeter-area relationship for the normal-phase 
clusters with fractal boundary. The chart is built for the trun- 
cated sampling containing 380 normal-phase clusters. The in- 
sets represent histograms of sampling for corresponding dis- 
tributions. 

accuracy of the statistical estimation, whereas the corre- 
lation coefficient (which becomes equal to 0.869) falls by 
six hundredth only. Therefore, we can conclude that the 
found estimate of the fractal dimension is robust. 



VI. THE PINNING GAIN FOR THE 
MAGNETIC FLUX TRAPPED IN FRACTAL 
CLUSTERS OF A NORMAL PHASE 

It has been revealed that the fractality of the cluster 
boundaries intensifies the pinning. As can be seen from 
Fig. [7J the decrease in the trapped magnetic flux at the 
same value of the transport current is less for larger frac- 
tal dimension. The relative change in the trapped flux 
A$/$, which, according to formula of Eq. J2J, can be 
calculated from Eq. (|18|) , defines the density of vortices 
n broken away from the pinning centers by the current i: 

B f B A$ 

n (0-^/(0 (») 

where B is the magnetic field, $o = he/ (2e) is the mag- 
netic flux quantum (h is Planck's constant, c is the ve- 
locity of light, e is the electron charge). The graph of 
the change in the trapped flux as a function of transport 
current, shown in Fig.0 coincides qualitatively with the 
curves obtained in experiments on the magnetization of 
YBCO films subjected to current pulsesi 68 i 69 Figure [7| 
also displays such an important property of supercon- 
ducting structure containing fractal clusters of a normal 
phase: the fractality intensifies the magnetic flux trap- 
ping, hindering its breaking away from pinning centers, 
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and thereby enhances the critical current which sample 
is capable to withstand, remaining in a superconducting 
state. The pinning amplification can be characterized by 
the pinning gain factor 

A$(L> = 1) 



201o gl 



dB (21) 



310 A$ (current value of D) 

which is equal to relative decrease in the fraction of mag- 
netic flux broken away from fractal clusters of fractal 
dimension D compared to the case of Euclidean ones. 
Dependencies of the pinning gain on the transport cur- 
rent for different fractal dimension at g — are shown 
in the inset of Fig. [7| The highest amplification (about 
10 dB) is reached when the cluster boundaries have the 
greatest possible fractality. Figure demonstrates that 
with increase in fractal dimension the trapped magnetic 
flux is changed less and less by the action of the transport 
current. The pinning gain of Eq. I|21|) characterizes the 
properties of a superconductor in the range of the trans- 
port currents corresponding to a resistive state (i > 1). 
At smaller current the total trapped flux remains nearly 
unchanged (see Fig.0) for lack of pinning centers of such 
small critical currents, so the breaking of the vortices 
away has not started yet. When the vortices start to leave 
the normal-phase clusters and move through the weak 
links, their motion induces an electric field, which, in 
turn, creates the voltage drop across the sample. There- 
fore, the passage of electric current is accompanied by the 
energy dissipation. As for any hard superconductor (that 
is to say, type-II, with pinning centers) this dissipation 
does not mean the destruction of phase coherence yet. 
Some dissipation always accompanies any motion of a 
magnetic flux that can happen in a hard superconductor 
even at low transport current. Therefore the critical cur- 
rent in such materials cannot be specified as the greatest 
non-dissipative current. Superconducting state collapses 
only when a growth of dissipation becomes avalanche- like 
as a result of thermo-magnetic instability. 

The principal reason of pinning enhancement due to 
the fractality of the normal-phase clusters lies in the 
fundamental properties of the critical current distribu- 
tion. Figure demonstrates the peculiarities of the frac- 
tal probability density specified by Eq. (JTSJ. As may be 
seen from this graph, the bell-shaped curve of the distri- 
bution broadens out, moving towards greater magnitudes 
of current as the fractal dimension increases. It means 
that more and more of the small clusters, which can trap 
the vortices best, are being involved in the game. Hence 
the number of vortices broken away from pinning centers 
by the Lorentz force is reducing, so the smaller part of a 
magnetic flux can flow. The shift of the critical current 
distribution towards higher currents can be described by 
dependencies of the mean critical current distribution on 
the fractal dimension, as it is shown in the inset of Fig. [S] 
The mean critical current obeys the strong super-linear 
law specified by Euler gamma function: 

(2+D)/2 , n s 

r(i--J (22) 
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FIG. 7: Effect of a transport current on the magnetic flux 
trapped in fractal clusters of a normal phase. The inset shows 
the pinning gain for the different fractal dimension of the 
cluster boundary. 



Figure |Hl clearly demonstrates that increasing the fractal 
dimension gives a growth of the contribution made by 
clusters of greater critical current to the overall distribu- 
tion, resulting just in enhancement of the magnetic flux 
trapping. 



VII. ANOMALOUS STATISTICAL 
PROPERTIES OF THE CRITICAL CURRENT 
FRACTAL DISTRIBUTION 

In the general case of gamma-distribution of the 
normal-phase cluster areas of Eq. (|14f> the probability 
density for the cluster areas has the form 



w(A) = 



A3 exp j-A/Ap) 

rfo + i; 



At 



(23) 



For further consideration it is convenient to introduce 
the dimensionless area of the cluster a, for which the 
distribution function of Eq. I|23|l can be rewritten as: 



w (a) 



(.9+1) 



3+1 



r G? + i) 



-a s exp(- (g+ l)a) (24) 



The mean dimensionless area of the cluster is equal to 
unity, whereas the variance is determined by one param- 
eter only: g\ — 1/ (g + 1). The critical current distribu- 
tion has the exponential- hyperbolic form of Eq. Ijl6|l , for 
which the mean critical current, instead of Eq. I)22[l. is 



g r(g + i-g/2) 
r(9 + i) 



(25) 
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FIG. 8: Influence of the fractal dimension of the boundary of 
the normal-phase clusters on the critical current distribution. 
The inset shows the dependencies of the mean critical current 
on the fractal dimension. 



The critical current distribution has the skew bell-shaped 
form with inherent broad "tail" extended over the region 
of high currents (see Fig. |HJ . The mean critical current of 
Eq. P5|l. just as mean current of Eq. in the special 
case of exponential distribution of the normal-phase clus- 
ters areas, is divergent in the range of fractal dimensions 
D> 2( 5 + l). 

The gamma distribution of the critical currents of 
Eq. HI 61 spreads out with some shift to the right as g- 
parameter decreases. This broadening can be character- 
ized by the standard deviation of critical currents 



1-D) (T{g+1- 



r (<? + i) 



r( ff - 



g/g) V 
1) J 



(26) 



The standard deviation grows nonlinearly with increase 
in the fractal dimension. The peculiarity of the distribu- 
tion of Eq. HI ftp is that its variance becomes infinite in 
the range of fractal dimensions D ^ g + 1. At the same 
time, the mode of the distribution modef(i) — (G/9) D / 2 
remains finite for all possible values of the fractal di- 
mension 1 ^ D ^ 2. The distributions with divergent 
variance are known in probability theory - the classic ex- 
ample of that kind is Cauchy distribution^ However, 
such an anomalous feature of exponential-hyperbolic dis- 
tribution of Eq. H16J) is of special interest, inasmuch as the 
current-carrying capability of a superconductor would be 
expected to increase in the region of giant variance. Then 
the statistical distribution of critical currents has a very 
elongated "tail" containing the contributions from the 
clusters of the highest depinning currents. 

The reason for divergence of statistical moments of 



the exponential-hyperbolic critical current distribution 
of Eq. H16II consists in the behavior of the cluster area 
distribution of Eq. H24|I .The form of this function takes 
essentially different shapes depending on the sign of g- 
parameter - from the skew unimodal curve at g > to 
the monotonic curve with hyperbolic singularity at zero 
point at g < 0. In the borderline case of g — 0, which 
separates these different kinds of the functions, the dis- 
tribution obeys an exponential law. It is just for negative 
values of ^-parameter that the mean critical current di- 
verges. The contribution from the clusters of small area 
to the overall distribution grows at g < 0. Since the clus- 
ters of small size have the least number of weak links over 
a perimeter, they can best trap the magnetic flux. There- 
fore, an increase of the part of small clusters in the area 
distribution of Eq. H24|) leads to a growth of the contribu- 
tion with high depinning currents made by these clusters 
in the critical current distribution of Eq. (|l(i|> . Just as a 
result of this feature the mean critical current diverges 
at g < 0. Nevertheless, the total area between the curve 
of the probability density and the abscissa remains finite 
by virtue of the normalization requirement. 

Obviously, the proper critical current cannot be in- 
finitely high as well as the clusters of infinitesimal area 
do not really exist. There is the minimum value of the 
normal-phase cluster area A mi which is limited by the 
processes of the film growth. So in YBCO based com- 
posites prepared by magnetron sputtering^ the sample 
value of minimum area of the normal-phase cluster has 
been equal to A m = 2070 nm 2 at mean cluster area A = 
76500 nm 2 , that corresponds to the lower bound of the 
dimensionless area of the cluster a m = A m j A — 0.027. 
In view of this limitation, we will describe the distribu- 
tion of the normal-phase cluster areas by the truncated 
version of the probability density of Eq. H24fl : 



w(a \a ^ a m ) 



h (a 



1- W(a m ) 



w (a) 



(27) 



where h (x) = 
function, and 

W (a m ) = 



for x ^ 
for x < 



is the Heaviside step 



/ \ j 7(.9+l,(3 + l)am) , osA 
w(a)da= ; r (28) 

r(s + i) 



is the truncation degree, which is equal to the probabil- 
ity Pr { Vo < a m } to find the cluster of area smaller than 
the least possible value a m in the initial population. In 
the expression of Eq. H28|) 7 (v, z) is the complementary 
gamma function. 

The expression of Eq. H27|l gives the conditional distri- 
bution of probability, for which all the events of finding 
a cluster of area less than a m are excluded. Thus the 
truncation provides a natural way to fulfil our initial as- 
sumption that the cluster size has to be greater than the 
coherence and penetration lengths. New distribution of 
cluster areas gives rise to the truncated distribution of 
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the critical currents: 

where z m = (G/ (g + 1) a m ) D is the upper bound of 
the depinning current, which corresponds to the cluster 
of the least possible area CLm . 

Then, instead of Eqs. (|26(l and JUJ, the standard de- 
viation and mean critical current are 

a * = G ° / r( g + l-£>,(g + l)a m ) 
V T(g + l,(g + l)a m ) 

/ r( g + l-£>/2,( g + l)a m ) y 

^ r( g + i,( g + i)a m ) ; 1 ' 



r = G oT( 9 + l-D/2,(g + l)a m ) 

T(g + l,(g + l)a m ) V 1 

For the truncated distribution of Eq. the possible 
values of depinning currents are bounded from above by 
the quantity i m , therefore the mean critical current as 
well as the variance do not diverge any more. Both of 
these characteristics are finite for any fractal dimensions, 
including the case of maximum fractality (D = 2). The 
truncation degree is related to the least possible area of 
the cluster by the equation (|28|l . 

The degree of truncation gives the probability measure 
of the number of normal-phase clusters that have the area 
smaller than the lower bound a m in the initial distribu- 
tion. If the magnitude of W(a m ) is sufficiently small (no 
more than several percent), then the procedure of trun- 
cation scarcely affects the initial shape of the distribu- 
tion and still enables the contribution from the clusters 
of zero area (therefore, of infinitely high depinning cur- 
rent) to be eliminated. It is interesting to note that the 
very similar situation occurs in analyzing the statistics 
of the areas of fractal islands in the oceanS The island 
areas obey the Pareto distribution, which also has the 
hyperbolic singularity at zero point that causes certain 
of its moments to diverge. For exponential distribution of 
the cluster areas, which is valid in the above-mentioned 
case of YBCO films^ the probability to find the clus- 
ter of area smaller than am in the sampling is equal to 
Pr{Va < a m } = 2.7% only. In principle, the truncation 
procedure could be made here, too, but there is no need 
for that, because the contribution of infinitesimally small 
clusters is finite at g — (and equal to zero at g > 0). 

As the fractal dimension increases, the critical current 
distribution broadens out (the variance grows), moving 
towards higher currents (the mean critical current grows, 
too). This trend is further enhanced with decreasing g- 
parameter. The most current-carrying capability of a 
superconductor should be achieved when the clusters of 
small size, which have the highest currents of depinning, 
contribute maximally to the overall distribution of the 
critical currents. Such a situation takes place just in the 



region of giant variance of critical currents. So far, the 
least magnitude of g -parameter (equal to zero) has been 
realized in YBCO composites containing normal-phase 
clusters of fractal dimension D = 1.44.— The critical 
currents of superconducting films with such clusters are 
higher than usual. 8,68,69 It would thus be expected to 
further improve the current-carrying capability in super- 
conductors containing normal-phase clusters, which will 
be characterized by area distribution with negative mag- 
nitudes of g -parameter, especially at the most values of 
fractal dimensions. 



VIII. CURRENT-VOLTAGE 
CHARACTERISTICS AND VORTEX DYNAMICS 
NEAR THE RESISTIVE TRANSITION 

By virtue of the fact that any motion of the magnetic 
flux causes the energy dissipation in superconductors, the 
question of how such a process could be prevented, or 
only suppressed, is of prime practical importance. The 
study of resistive state peculiarities leads to the conclu- 
sion that the cluster fractality exerts influence on the 
electric field induced by the flux motion. The found dis- 
tribution of the critical currents allows us to find the elec- 
tric field arising from the magnetic flux motion after the 
vortices have been broken away from the pinning centers. 
Inasmuch as each normal-phase cluster contributes to the 
total critical current distribution, the voltage V across a 
superconductor can be represented as the response to the 
sum of effects made by the contribution from each clus- 
ter. Such a response can be expressed as a convolution 
integral 

i 

— = (i-i')f(i')di' (32) 

Tf J 



where rf is the dimensionless flux flow resistance. The 
similar approach is used universally in all the cases 
where the distribution of the depinning currents takes 
place iSZ'ii'Z^ The subsequent consideration will be essen- 
tially concentrated on the consequences of the fractal na- 
ture of the normal-phase clusters specified by the distri- 
bution of Eqs. (|T5)) . l(T%|) . so all the problems related to 
possible dependence of the flux flow resistance rf on a 
transport current will not be taken up here. The voltage 
across a superconductor V and flux flow resistance Rf 
are related to the corresponding dimensionless quantities 
v and rf by the relationship: V/Rf = I c {v/rf). 

After substitution of the critical current distribution 
of Eq. I|16|l into the expression of Eq. I|32|) , upon integra- 
tion, we get the final expression for the current-voltage 
( V-I) characteristics in the general case of the gamma- 
distribution of the cluster areas: 
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In the special case of exponential cluster area distribu- 
tion (g=0) the general formula of Eq. i|33|) can be simpli- 
fied: 

^ = i exp (-Ci- 2 ' D ) - C D ' 2 T (\ - ^,Ci~^ D ) 

(34) 

where C = ((2 + D) /2) 2/D+1 . 

In extreme cases of Euclidean clusters and clusters of 
the most fractality the expression of Eq. 1341) can be fur- 
ther transformed: 

(i) Clusters of Euclidean boundary (D — 1 at g = 0): 

^zexp(^)-V3^erfc(^) (35) 

where erfc(z) is the complementary error function. 

(ii) Clusters of boundary with the maximum fractality 
(D = 2 at g= 0): 

~ = * cxp (~~) +4Ei ( _ ~) ( 36 ) 

where Ei (z) is the exponential integral function. 

The V-I characteristics calculated using the formulas 
of Eqs. 1231) - ESJ are presented in Fig. Two lines 
drawn for extreme case of D = 1 and D — 2 bound the 
region the V-I characteristics can fall within for any pos- 
sible values of fractal dimension. Figure El shows that 
the fractality reduces appreciably an electric field arising 
from the magnetic flux motion. This effect is especially 
strong in this range of the currents (1 < i < 3), where 
the pinning enhancement also has a maximum (see also 
Fig. UJ. Both these effects have the same nature, inas- 
much as their reason lies in the peculiarities of the critical 
current distribution. The influence of the fractal dimen- 
sion of the cluster boundary on the critical current distri- 
bution is demonstrated in Fig.[S] An increase of fractality 
causes a significant broadening of the tail of the distribu- 
tion / = f(i). Therefore, more of small clusters, which 
have higher critical currents of depinning, are being in- 
volved in the process. Hence the smaller number of vor- 
tices can move, creating the smaller electric field. In turn, 
the smaller the electric field is, the smaller is the energy 
dissipated when the transport current passes through the 
sample. Therefore, the decrease in heat-evolution, which 
could cause transition of a superconductor into a nor- 
mal state, means that the current-carrying capability of 
a superconductor containing such fractal clusters is en- 
hanced. 




1 2.3 4 5 
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FIG. 9: Current- voltage characteristics of superconductors 
with fractal clusters of a normal phase. 



Though the fractality of the clusters reduces the volt- 
age arising from the vortex motion in the range of the 
currents i > 1, the situation is quite different in the 
neighborhood of the resistive transition below the critical 
current. When i < 1, the higher the fractal dimension 
of the normal-phase cluster is, the larger is the voltage 
across a sample and the more stretched is the region of 
initial dissipation in V-I characteristic. The critical cur- 
rent i c is preceded by some initial region of the finite 
voltage drop starting with i on , so the resistive transition 
is not absolutely abrupt. The existence of this initial 
section on the V-I characteristic arises from the pecu- 
liarities of the fractal distribution in the range of small 
currents. It has been just a similar initial region of frac- 
tal dissipation has been observed in high-resolution mea- 
surements of dynamical resistance of HTS-normal metal 
composites iii 

The significant difference in V-I characteristic behavior 
below and above the resistive transition is related to de- 
pendence of free vortex density on the fractal dimension 
for various transport currents. The resistive characteris- 
tics provide additional information about the nature of a 
vortex state in superconductor with fractal clusters of a 
normal phase. Since the V-I characteristics of Eqs. 1)33(1 - 
(|34|) are non-linear, dc (static) resistance is not constant 
and depends on the transport current. The more conve- 
nient parameter is the differential resistance = dv/di, 
a small-signal ac parameter that gives the slope of the V- 
I characteristic. The corresponding dimensional quantity 
Rd can be found using the formulas Rd = TdRf/rf. Fig- 
ure EI| shows the graphs of the differential resistance as 
a function of transport current for superconductor with 
fractal normal-phase clusters. The curves drawn for the 
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FIG. 10: Dependence of the dc resistance on the transport 
current for superconductor with fractal normal-phase clus- 
ters. The dashed horizontal line r = 77 at the upper right 
corresponds to the viscous flow of a magnetic flux. 
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FIG. 11: Dependence of the free vortex density on the fractal 
dimension of the normal-phase clusters for different values of 
a transport current. The inset shows the free vortex density 
versus current for two values of a fractal dimensions. 



Euclidean clusters (D = 1) and for the clusters of the 
most fractal boundaries (D = 2) bound a region contain- 
ing all the resistive characteristics for an arbitrary fractal 
dimension. As an example, the dashed curve shows the 
case of the fractal dimension D = 1.5. The dependencies 
of resistance on the current shown in Fig. llUl are typical of 
the vortex glass, inasmuch as the curves plotted on a dou- 
ble logarithmic scale are convex and the resistance tends 
to zero as the transport current decreases, rd{i — > 0) — > 0, 
which is related to the flux creep suppressioniSSiSi A vor- 
tex glass represents an ordered system of vortices without 
any long-range ordering. At the same time, the vortex 
configuration is stable in time and can be characterized 
by the order parameter of the glassy stateiSiZi In the 
H- T phase diagram, mixed state of the vortex glass type 
exists in the region below the irreversibility line. The 
dashed horizontal line at the upper right of Fig. cor- 
responds to a viscous flux flow regime (r^ = rf =const), 
which can only be approached asymptotically. 

The resistance of superconductor is determined by the 
density n of free vortices broken away from pinning cen- 
ters of Eq. I|2U|) . In the special case of exponential cluster 
area distribution (g — 0) the free vortex density is equal 
to 

n= JL exp (-Ci- 2 / D ) (37) 

The more vortices are free to move, the stronger the in- 
duced electric field, and therefore, the higher is the volt- 
age across a sample at the same transport current. If we 
differentiate the V-I characteristic of Eq. I|34|l and com- 
pare with formula of Eq. i|37[l . we will see that the dif- 
ferential resistance is proportional to the density of free 
vortices: = (ry$ / B)n. Resistance of a superconduc- 
tor in the resistive state is determined by the motion just 



of these vortices. 

Figurc^Jdemonstrates dependence of the relative den- 
sity of free vortices n(D)/n(D — 1) (relatively to the 
value for clusters with Euclidean boundary) on the frac- 
tal dimension for different magnitudes of transport cur- 
rents. The vortices are broken away from pinning cen- 
ters mostly when i > 1, that is to say, above the resis- 
tive transition. Here the free vortex density decreases 
with increasing the fractal dimension. Such a behavior 
can be explained by the peculiarity of the critical cur- 
rent distribution of Eq. I|lt)l) . The curve of this distribu- 
tion broadens out, moving towards greater magnitudes 
of current as the fractal dimension increases. It means 
that there are more clusters of high depinning current in 
superconductor. The smaller part of the vortices is free 
to move, the smaller the induced electric field. The fact 
that the fractality of cluster boundary enhances pinning 
is also demonstrated in Fig. ^3 where the resistance de- 
creases with increasing fractal dimension above the resis- 
tive transition. The relative change in free vortex density 
depends on the transport current (see inset in Fig. Ill|) 
and, in the limiting case of the most fractal boundary 
D = 2 reaches a minimum for i = 1.6875 (this curve 
should goes below others in Fig. Illfl . That corresponds 
to the maximum pinning gain and the minimum level of 
dissipation. 

In the range of transport currents below the resistive 
transition (i < 1) the resistance as well as the free vortex 
density increase for the clusters of greater fractal dimen- 
sion (see Figs. I1UI and Illfl . Such a behavior is related to 
the fact that the critical current distribution of Eq. 1|16|) 
broadens out covering both high and small currents as the 
fractal dimension increases. For this reason, the break- 
ing of the vortices away under the action of transport 
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current begins earlier for the clusters of greater fractal 
dimension. In spite of sharp increase in relative density 
of free vortices (Fig.HT|l. the absolute value of vortex den- 
sity in the range of currents involved is very small (much 
smaller than above the resistive transition). So the vortex 
motion does not lead to the destruction of the supercon- 
ducting state yet, and the resistance remains very low. 
The low density of vortices at small currents is related to 
the peculiarity of exponential-hyperbolic distribution of 
Eq. mJJ|. This function is so "flat" in the vicinity of the 
coordinate origin that all its derivatives are equal to zero 
at the point of i = 0: d k f(0)/di k = for any value of k. 
This mathematical feature has a clear physical meaning: 
so small a transport current does not significantly af- 
fect the trapped magnetic flux because there are scarcely 
any pinning centers of such small critical currents in the 
overall statistical distribution, so that nearly all the vor- 
tices are still pinned. This interval corresponds to the 
so-called initial fractal dissipation regime, which was ob- 
served in BPSCCO samples with silver inclusions as well 
as in polycrystalline YBCO and GdBCO samples. 11 



the fractal properties of the clusters and dynamics of the 
trapped magnetic flux is established. The fractal distri- 
bution of the critical current is obtained. It is found that 
the fractality of cluster boundary strengthens the flux 
pinning and thereby hinders the destruction of supercon- 
ductivity by the transport current. V-I characteristics 
of fractal superconducting structures in a resistive state 
are obtained. It is revealed that the fractality of the 
boundaries of the normal-phase clusters reduces the elec- 
tric field arising from magnetic flux motion, and thereby 
raises the critical current of a superconductor. The frac- 
tal properties of the normal-phase clusters significantly 
affect the resistive transition. This phenomenon is re- 
lated to the features of the fractal critical current dis- 
tribution. The resistance characteristics correspond to a 
mixed state of the vortex glass type. An important re- 
sult is that the fractality of the normal-phase clusters en- 
hances pinning, thus decreasing the resistance above the 
resistive transition and increasing the current-carrying 
capability of a superconductor. 



IX. CONCLUSION 
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